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ABSTRACT 

In a universe where AGN feedback regulates star formation in massive galaxies, 
a strong correlation between these two quantities is expected. If the gas causing star 
formation is also responsible for feeding the central black hole, then a positive corre¬ 
lation is expected. If powerful AGNs are responsible for the star formation quenching, 
then a negative correlation is expected. Observations so far have mainly found a mild 
correlation or no correlation at all (i.e. a flat relation between star formation rate 
(SFR) and AGN luminosity), raising questions about the whole paradigm of “AGN 
feedback”. In this paper, we report the predictions of the GALFORM semi-analytical 
model, which has a very strong coupling between AGN activity and quenching of star 
formation. The predicted SFR-AGN luminosity correlation appears negative in the 
low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive 
in the regime of the brightest AGN. Our predictions reproduce reasonably well recent 
observations by [Rosario et al-l yet there is some discrepancy in the normalisation of 
the correlation at low luminosities and high redshifts. Though this regime could be 
strongly influenced by observational biases, we argue that the disagreement could be 
ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the 
galaxies that dominate the regime where the observations imply a weak correlation 
are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these 
galaxies retain high enough molecular hydrogen contents to maintain relatively high 
star formation rates and strong infrared emission. 
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1 INTRODUCTION 

A clear understanding of how galaxies transition from star¬ 
forming disks to passive spheroidals is one of the open prob¬ 
lems in the current paradigm of galaxy formation. One of 
the most widely accepted theories requires a large energy 
injection into the cores of massive galaxies; such energy will 
then heat up the cold gas in the center and prevent the ac¬ 
cretion of new gas onto the galaxy. The most likely source 
for such energy are active galactic nuclei (AGN), powered by 
matter accretion onto super massive black holes (SMBHs) 
at the centers of galaxies ( Silk fc Rees|1998| ). 

In the last decade both semi-analytical models and nu¬ 
merical simulations have tried to incorporate such an ef¬ 
fect into the modelling of galaxy formation. Semi-analytical 
models have showed that if the energy from AGN is coupled 
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to the hot and cold gas, it is indeed possible to halt star for- 


mation and create passive galaxies 

Kauffmann & Haehnelt 
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Croton et al. 2006| Somerville et al. 2008| Fanidakis et al.j 
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I. Numerical simulations have also 


helped to shed light onto the interplay between a galaxy 


and its SMBH (e.g. Springel et al. |2005a 

Di Matteo et al. 

20051 

Robertson et al. 2006 

ISijacki et a 

. |2007| Hopkins 

et al.| 

2007 Di Matteo et a 

.| 20081 lOkamoto et al. 2008| 

Booth & Schaye||2009 More recently, large fully cosmolog- 


ical simulations including the effect of AGN feedback have 
been quite successful in reproducing properties of observed 
galaxies, including the luminosity function and star forma¬ 
tion rate (Vogelsberger et al.||2014 and Schaye et al.||20l5|. 


There is general agreement that AGN feedback is a needed 
ingredient in galaxy formation and evolution. 

These studies imply a strong correlation between the 
evolution of a galaxy and its central black hole. A hint to 
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such a correlation might reside in the so called M-sigma rela¬ 
tion, namely the observed correlation between the velocity 
dispersion of the bulge of the galsixy and the mass of the 


supermassive black hole hosted by it (Ferrarese & Merritt 
|2000[ IGebhardt et al.||2000[ but see mke fc Maccio||2011 
for a different explanation). 

Another possible piece of evidence of a causal link be¬ 
tween the rising of AGN activity and the quenching of a 
galaxy can be obtained by studying the observed correla¬ 
tion between star formation rate (SFR) and AGN luminos¬ 


2012 ). 


ity (Shao et al. 2010 Lutz et al. 2010 and Harrison et al. 


Recently, [Rosario et al.|[2012| used deep infrared and 
X-ray observations in the COSMOS, GOODS-NORTH and 
GOODS-SOUTH fields to estimate the SFR and AGN lumi¬ 
nosity, respectively. They found a weak correlation between 
SFR and AGN activity at all redshifts, suggesting little con¬ 
nection between SF and BH growth in these systems, espe¬ 
cially at high redshift and low luminosity. [Mullaney et al.| 
|2012b| and [Stanley et al.|2015|find a similar absence of cor¬ 


relation (but see Mullaney et al.|2012a and Rodighiero et al. 


|2015| for a possible correlation of SMBH growth rate and 
SFR). 

This quite unexpected lack of correlation between SFR 
and AGN luminosity has been used to put constraints on 


the triggering of AGN activity and black hole growth | Neis- 


tein & Netzer 2014), possibly suggesting strong variability 


in AGN activity on time scales shorter than those typical 


of star formation, which are of the order 100 Myrs (Hickox 
|et al.|2014| ). 

In this paper, we take a deeper look into the expected 
correlation between SFR and AGN activity. We present pre¬ 
dictions from the GALFORM semi-analytical model of galaxy 
formation ( Cole et al.||200d l in which AGN feedback is the 
key ingredient to explain the fading of star formation in 
massive galaxies. GALFORM has been successful in reproduc¬ 
ing a number of fundamental relations of galaxy evolution 
and structure: the luminosity and stellar mass functions 
of galaxies, the AGN diversity and evolution, the evolu¬ 
tion of Lyman-Q emitters (LAEs) and Lyman-break galax¬ 
ies (LBGs), the number counts of sub-millimeter galaxies 
(SMGs), as well as the HI and HII mass functions ([Baugh 


et al.j[2005[Bower et al.[[20(I6[ [Fanidakis et al.[[2011[ [Lac^ 

et al.j20'rT and Lagos et al.|2012 I. 


Using the version of GALFORM presented in [ Fanidakls[ 
[et al.[ [2012[ and assuming full co-evolution between the 
galaxy and its central BH, we study the predicted correlation 
between SF and AGN luminosity and compare it to obser¬ 
vations and previous models of this relation. We stress that 
this correlation is not used as a constraint on the model, but 
presents a pure prediction of a model constrained by other 
observations. 

The paper is organised as follows. In Section 1, we in¬ 
troduce the GALFORM semi-analytical model and describe the 
most important model ingredients that are relevant for this 
study. In Section 2, we present the predictions of the model 
for the FIR luminosity as a function of AGN activity and 
compare to the observations. In Section 3 we look at the 
molecular gas content of our galaxies to understand the ma¬ 
terial that is causing star formation. Finally, in Section 4, 
we discuss our findings and present our conclusions. 


2 THE MODEL 


Semi-analytic modelling in galaxy formation is used to simu¬ 
late large numbers of galaxies in a computationally efficient 
manner, thus enabling statistical insight into the predictions 
of different galaxy formation models. It uses our understand¬ 
ing of cosmological structure formation in N-body simula¬ 
tions and adds a set of coupled differential equations to de¬ 
scribe the physical processes of galaxy formation and evolu¬ 
tion, making possible the analysis of a wide variety of galaxy 
and BH properties. GALFORM, one of the most extensively ap¬ 
plied semi-analytical models, has been improved upon over 
time ( Gole et al.|2000 Baugh et al.|2005 Bower et al.|2006 


Fanidakis et al.|2011 Lagos et al.|2011 1 to include processes 

and adjustments that describe more accurately the forma¬ 
tion and evolution of galaxies. Among the most important 
physical processes that are modelled in GALFORM are: i) the 
formation and evolution of DM haloes in the AGDM cos¬ 
mology, ii) gas cooling and disk formation in DM haloes, 
iii) star formation, supernovae feedback and chemical en¬ 
richment, iv) BH growth and AGN feedback, v) and the 
formation of bulges through galactic disk instabilities and 
galaxy mergers. 

The model is successful in reproducing many observa¬ 
tions, including the luminosity and stellar mass functions 
of galaxies, the number counts of submillimeter galaxies, 
the evolution of Lyo and Ly-break galaxies, H and H 2 mass 


functions and the AGN diversity and evolution (Baugh et al. 

2005[ Bower et al.|2006 Orsi et al. 

2008[ [Kim et al. 2011[ 

Lacey et al.| 2011[ [Lagos et al.|2011 

Lagos et al.||2012[ [La- 

gos et al.|[2014a[ [Lagos et al.||2014b 

Fanidakis et al. 2011[ 

h’anidakis et al.|2012[ h’anidakis et al 

2013a[ h’anidakis et al. 

2013b [Gonzalez-Perez et al. 2013[ 

Gonza ez-Perez et al. 

2014). The model provides a very good fit to the HER- 


SGHEL far-IR galaxy luminosity functions (LFs) and num¬ 
ber counts (Lacey et al. in prep.; see also Lacey et al.|201(3 ) 
and X-ray LFs of AGN in the 2 = 0 — 3 universe ( [Fanidakis 
et al. 2012[ Fanidakis et al. in prep.) and therefore it is 
an ideal tool for studying the correlation between SFR and 
AGN activity. 


For the purposes of this analysis, we use a version 
of the [Fanidakis et ah] ( [2012| model, updated to the cos¬ 
mological parameters estimated by the 7-year data release 
of WMAP ( [Komatsu et al. 2011 and the Lacey et al. in 
prep. GALFORM model). The fundamental predictions of the 
model for the AGN properties are described in [Fanidakls| 
et al. ( [2013a|b[ see also Fanidakis et al. in prep.). The 


predictions presented in this analysis are calculated using 
merger trees extracted from the DM only N-body simula¬ 
tion Millennium WMAP-7 (Lacey et al. in prep.). The Mil¬ 
lennium WMAP-7 simulation has almost the same mass res¬ 
olution (8.61 X 10® M 0 ), particle number (10^) and box size 


(500Mpc) as the Millennium simulation ( Springel et al. 
[2005b| and differs only in the background cosmology (which 
is in agreement with WMAP7 results). We now describe in 
the rest of this section the main model processes and in¬ 
gredients that are essential for understanding the GALFORM 
predictions for the SFR and AGN activity. 
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SFR- AGN luminosity correlation 3 


2.1 The two regimes of black holes growth 

GALFORM assumes two regimes of black hole (BH) growth 
in galaxies. The first one is the “starburst” regime, where 
a galaxy experiences a starburst and AGN activity with 
high star formation and BH accretion rates. The high ef¬ 
ficiency both in SF and BH growth of this regime is due to 
efficient gas cooling in DM haloes with masses lower than 
5 X 10^^ Mq. The onset of the starburst mode is triggered by 
either a merger or by a disk instability. In this regime, a frac¬ 
tion of the cold gas that turns into stars is accreted directly 
onto the BH. When only this regime is active, we expect a 
positive correlation between star formation rate (SFR) and 
active galactic nucleus (AGN) luminosity. 

In the “hot-halo” regime, the growth of BHs is tightly 
linked to the AGN feedback mechanism and the suppres¬ 
sion of gas cooling in haloes typically more massive than 
5 X 10^^ Mq. The gas feeding the BH during this mode is 
assumed to originate directly from the hot halo around the 
galaxy. The resulting accretion power is then coupled via a 
jet to the thermodynamical properties of the gas in the host 
halo and suppresses cooling, if the available heating power 
exerts the cooling luminosity of the gas, I/cooi. Galaxies that 
are subject to AGN feedback still exhibit some SF, due to 
remaining cold gas and new gas brought in by mergers, but 
it is expected to be lower than in actively SF galaxies. In 
this regime, we expect a negative correlation between SFR 
and (AGN) luminosity. 

In addition to the starburst and hot-halo modes of BH 
growth, BH-BH mergers during galaxy encounters also con¬ 
tribute to the growth of BHs. However, this growth mode 
only redistributes the BH mass and does not add new 
baryons to the BHs. 


2.2 Bolometric accretion luminosity 

The gas accreted during the starburst mode is converted 
into an accretion rate, M, by assuming that the accretion 
duration is proportional to the dynamical timescale of the 
host spheroid, 

M=^^. (1) 

J qT'bulge 

Here M^cc is fixed for every galaxy to 0.5 percent of the mass 
that turns into stars during a starburst, tdyn is the dynamical 
timescale of the host spheroid and /q is a proportionality 
factor set to 10 in |Fanidakis et ar] ( |2012[ ). In contrast, in the 
hot-halo mode, the accretion rate onto the BH is calculated 
directly from the cooling properties of the host DM halo, i.e. 

M=^, (2) 

CkinC 

where tkin is the average kinetic efficiency of the jet dur¬ 
ing the AGN feedback. Lcooi is the quasi-hydrostatic cool¬ 
ing luminosity of the halo. This is chosen assuming that the 
flow will balance heating and cooling in the hot halo mode 
(i.e. if the Eddington ratio is sufficiently large.) This accre¬ 
tion mode is associated with early type galaxies in massive 
haloes, with relatively low SFRs. Since the accretion pro¬ 
cess is responsible for shutting down SF in the host galaxy, 
a negative correlation between AGN luminosity and SFR is 
expected in this mode. 


The bolometric luminosity of the accretion flow, Lboi, 
is calculated by assuming the Shakura-Sunyaev thin disk 
model I Shakura & Sunyaev |[l973t : 


Ebol.TD = eMc . 


( 3 ) 


We assume that this solution is valid for accretion rates 
higher than 1 percent of the Eddington accretion rate, 
i.e. m = M/M^dd ^ 0.01, where MEdd is defined 

as L-Edd/eF- For lower accretion rates, the advection- 
dominated accretion flow (ADAF) thick-disk solution is used 
( [Narayan fc Yi|[i994[ |Mahadevan||1997| , 

Eboi,ADAF = 0.44 eMc^. (4) 

When the accretion becomes substantially super-Eddington 
(Eboi ^ nLEdd), the bolometric luminosity is limited to 
( [Shakura fc Sunyaev||1973| : 

Eboi(^ vLEdd) = Till + ln(m/77)]LEdd, (5) 


where rj is an ad hoc parameter, which we choose equal to 
4, to allow a better modelling of the bright end of the AGN 
luminosity function (see Fanidakis et al. 20121. However, 
the accretion rate, rh, is not restricted if the flow becomes 
super-Eddington. 


2.3 Star formation rate 


Lagos et al. (20111 recently revised the original formula¬ 
tion of the SF law in GALFORM ( Cole et al.pOOO' l, in favour 
of the observationally motivated (and constrained) Blitz & 
Rosolowsky ( Blitz fc Rosolowsky||20(i6 1 empirical law. The 
Blitz & Rosolowsky law is motivated by UV, FIR and mil¬ 
limetre observations of SF and molecular gas in spiral galax¬ 
ies and assumes that the surface density of star formation 
rate (SFR) is proportional to the surface density of molecu¬ 
lar hydrogen in the disk. 


EsFR = TTSFEmol. 


( 6 ) 


The proportionality factor, vsf, is given as an inverse time 
scale and its value is closely constrained by observations. At 
every timestep, GALFORM calculates the molecular-to-atomic 
hydrogen ratio in the disk, Rmoi, which is expressed as a 
power law of the internal hydrostatic pressure of the disk 
( Elmegreen[]1993 1: 


^ E(H2) ^ /Pext\ 

S(H) [Po J 


( 7 ) 


Eq.jej is then re-written in terms of the cold gas surface 
density, Egas, as: 


EsFR — TTgp J'jjjQ^Egas, (8) 

where /moi = Pmoi/(l + Pmoi). Integration of Eqn. [^over 
the assumed exponential surface density profile of the gas 
gives the SFR of the galaxy. 

When the galaxy experiences a burst of SF, triggered 
either by a galaxy-galaxy merger or a disk instability, the 
model assumes that the entire cold gas reservoir of the 
galaxy (atomic and molecular) is converted into stars. The 
SF timescale in starbursts is finite and proportional to the 
dynamical timescale of the host bulge. At a given time, we 
estimate the total SFR in the galaxy as the sum of the qui¬ 
escent SFR and burst SFR. 
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2.4 Stellar and dust emission 


GALFORM includes a self-consistent model for calculating the 
stellar emission in every galaxy and its absorption and re¬ 
emission in the mid and far-IR (8 — 1000 p.m) and sub-mm 
wavelengths by dust. The model is similar in outline to the 
spectro-photometric model GRASIL ( Silva et al.||1998 l, with 
some simplihcations that speed-up signihcantly the calcula¬ 
tion. The stellar emission in the model is calculated based 
on the stellar population synthesis model of |Bressan et ah] 
119981 and assuming that the stars have an axisymmetric 
distribution in the disk and bulge. Dust in the galaxies is 
assumed to be a two-phase medium, i.e. diffuse low-density 
dust in the ISM and dense dust clouds enshrouding star 
forming regions. The amount of dust in each galaxy is de¬ 
termined by the total mass and metallicity of cold gas. The 
attenuation of stellar light by the dust is computed by inter¬ 
polating the tabulated radiative transfer models of |Ferra^ 
let al. 119991 The FIR emission from the dust is calculated as¬ 
suming that the cloud and diffuse dust components re-emit 
the absorbed starlight (i.e. the difference between the stel¬ 
lar and dust attenuated SED, integrated over wavelength) as 
a modified blackbody. This constrains the dust temperature 
(constant within a galaxy and different for each component), 
which is then used to calculate the SED of each dust com¬ 
ponent. The total dust SED is then found by summing over 
the SEDs of the two components. This technique works well 
for wavelengths greater than ~ 60pm. 


The composite SED of stellar and dust emission of 
galaxies in GALFORM provides a very good fit to galaxy 
number counts and luminosity function in the FIR (see 
Lacey et al. in prep.). Here we are interested in compar¬ 
ing the predictions of the model directly to far-IR obser¬ 
vations of AGN with the Herschel/PACS instrument, at 
a mean rest-frame wavelength of 60pm. For model galax¬ 
ies predicted by GALFORM, the monochromatic luminosity 
I'Leotim = i^Li/(60pm) is typically ~ O.OSLir, where Air is 
the total FIR luminosity, defined as the integrated luminos¬ 
ity between 8 and 1000 pm. This wavelength range includes 
the entire emission from dust, while excluding the intrinsic 
stellar emission. In the rest of the analysis, the GALFORM pre¬ 
dictions will be presented at 60pm, unless the comparison 
with observational data requires a different band or the total 
dust emission. We note that the galaxy SED predicted by 
GALFORM does not include any AGN emission. Therefore, the 
dust emission predicted by the model is purely reprocessed 
stellar emission. 


Finally, we note that the SFR is usually traced by emis¬ 
sion in the FIR as the SFR is not a direct observable. As 
a hrst test, we check how good of a proxy this emission 
is for the actual SFR in our model. Fig. shows the me¬ 
dian IR luminosity at 60pm (blue line) and its 10th and 
90th percentile (green shaded areas) as a function of SFR at 
z = 0 and 2.1. For high redshifts, the correlation is tighter, 
but even at a = 0 it is still present over several orders 
of magnitude. The positive dependance of Air on SFR is 
a consequence of the higher UV luminosity caused by the 
young stars. This in turn increases the dust heating and 
IR emission. At high SFRs, the 60pm emission seems to be 
directly proportional to the SFR. At intermediate and low 
SFRs (SFR < lO'^ — 10^ MQyr~^) the correlation deviates 
mildly from a slope of unity mainly due to a somewhat lower 




Figure 1. Predicted FIR emission at 60 pm as a function of 
star formation rate (SFR) at z = 0 and z = 2.1. The solid lines 
show median values, and the shaded areas the 10th and 90th 
percentiles. The red dashed line shows a line with a slope of unity 
and an arbitrary normalization. 


dust extinction or lower dust temperature. Overall, we can 
conclude that the FIR emission strongly couples to the SFR, 
thus, IR luminosity can be used as a proxy for SFR. We will 
henceforth use the FIR emission at zzAeo^jm as a proxy of 
star formation (SF). 


3 THE FIR-AGN LUMINOSITY 
CORRELATION 

In this study we are interested in comparing the predicted 
and observed FIR—AGN luminosity correlations. Our final 
aim is to check what correlation of observable quantities 
is predicted by a model where AGN feedback is the driv¬ 
ing mechanism for quenching star formation in high mass 
galaxies. 


3.1 Model Predictions 

We show our first predictions for the FIR—AGN luminosity 
correlation in Fig.[^ where we plot the distribution of galax¬ 
ies on the Aboi — i^Aeo^m plane. These galaxies represent a 
small subset (4, 000) of the total sample (usually of the order 
of 10® galaxies) and are randomly selected from the model 
output at 2 = 2.1. We split our sample into galaxies undergo¬ 
ing active starburst AGN activity (cold-gas accretion, blue 
points) and galaxies whose BHs are in the hot halo mode 
(hot-gas accretion, red points). We also divide the sample 
by stellar mass, with galaxies of M* > marked by 

a larger symbol, to aid the reader in finding the trends in 
galaxies that are more likely to be observed. Interestingly, 
low luminosity objects accreting in the hot-halo mode show 
a negative trend of i^Aeo/jm with AGN luminosity. In con¬ 
trast, high luminosity AGN, powered through the starburst 
mode, show a clear positive trend. The two different trends 
are a manifestation of the fact that each of the two accretion 
modes is linked to a different regime of SF efficiency. 

The different trends are due to the two different modes 
in which a BH accretes gas. Large AGN luminosities are 
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Figure 2. Predicted FIR emission at 60 |J.m, i/Leo^im, as a func¬ 
tion of bolometric AGN luminosity, Lboh for 4,000 randomly se¬ 
lected galaxies a,t z = 2.1 (both centrals and satellites). Galaxies 
in which the central BH experiences starburst-mode activity are 
shown with blue, while those with BHs growing during the hot- 
halo mode with red. The point size indicates stellar masses below 
(small) and above (big) 1 O^'^M 0 . 


due to accretion of cold gas and mostly happen in haloes 
with masses below ~ 5 x 10^^ Mq, as shown in 


Fanidakis 


et al. (2013b I. In this case, the triggering mechanism (galaxy 


mergers or disk instabilities) of cold-gas flows that feed the 
central BH is also responsible for a burst of SF in the host 
galaxy, hence the positive correlation. Low luminosity AGN 
mainly live in high mass haloes (> 10^® M©), where the gas 
fuelling the BH is accreted directly from the hot halo around 
the host galaxy (see Section]^. Accretion activity in this 
mode is tightly linked to the AGN feedback mechanism and, 
thus, is responsible for the suppression of gas cooling and 
SF. The quenching nature of this mode gives rise to a nega¬ 
tive correlation. In summary, the intrinsic uLeofim vs Lagn 
correlation predicted by GALFORM is entirely shaped by the 
physics of each of the accretion modes that are responsible 
for growing the central BH. 


3.2 Comparison with Data 



Figure 3. Top: Average FIR emission at 60p.m as a function of 
AGN bolometric luminosity, Lboli for all galaxies in our model. 
Different lines indicate different redshifts as indicated by the key. 
The shaded area represents the 10-90 percentiles for the model at 
z=2.I (pink), but percentiles are similar at all redshifts. Points 
with error bars show the observed — Tbol correlation in 

Rosario et al. (2012), where the z = 0 data are from Swift-BAT. 
The solid-dashed line represents the correlation line connecting 
various AGN datasets (SDSS type-II AGN with 0.1 Si z ^ 0.7 and 
Lbol Si 10"*^ erg s“^, and Spitzer/IRS type-I quasars at z ~ 0.1 
and z = 2 — 3) on the SFR — Lbol plane in Netzer (2009). Bot¬ 
tom: The predicted SFR-Lbol correlation as solid lines at different 
redshifts. We convert the Rosario et al. data points into SFRs ac¬ 
cording to Neistein & Netzer (2014). 


To compare our predictions with observational data, we 
show in Fig. the average izLeo/rm vs. Lagn relation now 
for the whole galaxy sample and calculated in five different 
redshifts bins (solid lines). For the calculation of the average 
I'Leo/jm we use the complete sample of galaxies at a given 
redshift without imposing any FIR luminosity cut. Our pre¬ 
dictions are compared with the observational results from 
the BAGS Evolutionary Probe (PEP) Herschel survey for the 
mean EIR emission of X-ray selected AGN at 0.4 ^ z ^ 2.1 
I Rosario et al. 2012| ). Also shown is the average FIR emission 
of local [z = 0) X-ray selected AGN from the Swift-BAT 
sample i Gusumano et ah] 20101. The straight dashed line 
has a slope of 0.8 [uLeofim oc Lagn) and shows the corre¬ 
lation line connecting various observational datasets on the 
SFR — I/AGN plane in |Netzer|r2009| These datasets include 
type-H SDSS AGN (0.1 ^ z ^ 0.7 and Lboi > lO'^^erg s“^) 
and Spitzer type-I quasars at z ~ 1 ( [Netzer et al.|[?007[ ) 
and z = 2 — 3 (Lutz et al.||2008|. Interestingly, the predic¬ 


tions of GALFORM for the izLeo/jm — Lagn correlation show the 
same behaviour at all redshifts: for AGN with bolometric lu¬ 
minosities below 10"^® ergs/s there is a negative correlation 
between SF and Lagn, while the trend is reversed (e.g. pos¬ 
itive correlation) for higher luminosities. 

Our model agrees well with the Swift-BAT data at 
z = 0. In addition, the predicted trend for high luminosity 
AGN, i.e. Lagn ^ 10“^® erg s“^, follows the slope of the cor¬ 
relation line by Netzer, although the model flattens mildly in 
the regime of the brightest quasars. Moreover, the GALFORM 
results at Lagn ^ erg s“^ follow the observational data 
from Rosario et al. . Both [Rosario et al.||2^012| and [Stanley[ 
[et al.[2015[note an increase in izLeoMm at higher Lagn, albeit 
the increase is weaker in [Stanley et al.|2015| Our model does 
predict a moderate increase of SFR towards higher Lagn- 
In our case, this is due to the fact that the high AGN lu- 
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model proposed by Hickox et al. (even though the SFR 
and AGN activity are fully and directly coupled in GALFORM 
without any AGN variability). The Lir — Tagn correlation 
predicted by the model is entirely shaped by the starburst 
mode, (this is in agreement with the observed correlation 
between SMBH growth rate and the SFR reported by |Mul-| 
laney et al.||201^|. In contrast to the predictions shown in 


ww we now find a monotonically increasing correlation. 
The average Lagn Inminosity is dominated by those galax¬ 
ies undergoing starburst-mode AGN activity, and thus the 
slope of the correlation reflects the slope of the starburst- 
mode branch in Fig.[^ This is because the typical AGN lu¬ 
minosities of these galaxies, and also their number density, 
are much higher than the ones of galaxies in the hot-halo 
mode branch. 


Figure 4. Average AGN bolometric luminosity, Lboi, as a func¬ 
tion of Ljr (8 — 1000|xm), to allow a direct comparison with 
the observational data. Different redshifts are shown by different 
colours, as indicated by the key. Points with error bars show the 
observational data by Rafferty et al. (2011), Chen et al. (2013), 
and Symeonidis et al. (2011). The dot-dashed line represents 
the predictions of the theoretical model presented in Hickox et 
al. (2014). 


minosity end is dominated by the starburst mode where a 
strong correlation between SFR and Lagn is assumed (see 
section EB- 

However, at lower luminosities and higher redshifts the 
model diverges from observations. There is a pronounced 
dip in the model at around Tagn ~ 10"^® erg s“^, apparently 
where the model’s two modes coincide (see figure]^, poten¬ 
tially an artifact of the abrupt transition between the two 
modes. The negative trend at low luminosities is not visible 
at all in the data from Rosario et al. (2012), which are prac¬ 
tically flat at all redshifts. This observed fiat correlation was 
used in previous studies (Rosario et al.| Mullaney et al. I as 


a hint of a possible disjoint evolution of BHs and their host 
galaxies, contrary to what is commonly assumed in many 
models of galaxy formation. Here we find that in GALFORM, 
this lower luminosity regime is purely shaped by the strong 
correlation between AGN activity and galaxy evolution. 

We also show in the same plot (lower panel) the correla¬ 
tion between AGN luminosity and average SFR for the same 
sample of galaxies shown in the top panel. The Herschel ob¬ 
servational data have now been divided by 1.9 x 10“^^ erg s“^ 
to convert the observed the fiOpm flux into a rough esti¬ 
mation for the SFR (Neistein & Netzer 20141. The model 
predictions show a behaviour similar to the vLeo/im — Lagn 
correlation shown in the top panel, though the disagreement 
with the observational data is now more evident. We note 
however, that our estimate of the observed SFR is only a 
crude approximation, and therefore this comparison with 
the model predictions is only for illustrative purposes. 

Finally, we plot in Fig. |^the average AGN luminosity 
Tagn as a function of the total FIR luminosity Lir. Our 


predictions are compared to recent observational data (Raf- 


ferty et al.||20lf| |Symeonidis et ^|2011| and [Chen et al. 


2013) and the AGN variability model proposed by 


Hickox 


et al.|[^14[ Our model predictions are in moderate agree¬ 


ment with the observational results, as well as with the 


3.3 Properties of AGN hosts 

The Lagn — Lib. correlation can in principle provide in¬ 
sights into the host properties of AGN and thus, impose 
constraints on theoretical models of galaxy formation. For 
example [Rosario et al.| suggest that the flatness of the cor¬ 
relation, LAGN-vLeofim, at low AGN luminosities and its 
steep evolution at higher luminosities indicates two regimes 
of AGN activity. What our results show is that these two 
AGN regimes happen in very different galaxy populations. 

AGN at the bright end of the Lagn — I'Leo/j.m corre¬ 
lation shown in Fig. are found in bursty systems expe¬ 
riencing intense SF. The hosts of these luminous AGN are 
gas rich systems living in ~ DM halos, and are 

similar to actively star forming galaxies. In contrast, in the 
faint Lagn region of the plane we find galaxies whose cen¬ 
tral BHs accrete predominately via the hot-halo mode. As 
already mentioned, this mode is particularly prominent in 
more massive systems (> 10^^ Mq), where AGN feedback 
efficiently suppresses gas cooling and SF. The typical host 
of an AGN in that region would be a spheroidal system that 
at low redshifts resembles an early-type galaxy. 

The stellar masses of the galaxies populating the 
Lagn — Lir plane span a wide range of values (see Fanidakis 
et al. 2013b). We show this in Fig.j^ where we plot the SFR 
as a function of stellar mass at 2 = 0 for active and inac¬ 
tive galaxies in our model. We define as inactive all galaxies 
with Lagn < 10"^^ erg/s. At the same time, we split active 
galaxies into faint (10"^^ < Lagn < 10"*^ erg/s) and bright 
AGN (Lagn > 10’’^ erg/s). The transition luminosity from 
faint to bright AGN approximately marks the break in the 
Lagn — vLqoij.ui correlation shown in Fig. We refer the 
reader to [Lagos et al.j ( |2011| ) for the overall properties of 
galaxies on the Mgtar — SFR plane in GALFORM. 

What is immediately evident from Fig. [^is that AGN 
and inactive galaxies populate the same regions on the 
Mstar — SFR plane. Bright AGN are predominately found 
on the main sequence, though the scatter is strong. AGN in 
this sample are powered by accretion during the starburst 
mode. On the other hand, faint AGN are found both on the 
main sequence, but also immediately below it, i.e. in the 
region of passive galaxies. AGN hosts in the passive regime 
are typically quenched and have low SFRs. In the faint AGN 
sample we find AGN powered by both the starburst (AGN 
in the main sequence) and the hot-halo mode (AGN in the 
passive regime). Overall, we find no particular feature, e.g. 
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Figure 5. Stellar mass - star formation rate plane at z = 0 
for an unbiased subset of model galaxies in three bins: inactive: 
I'AGN < 10'^^ erg/s, faint AGN: 10"^^ < I/agn < 10^^ erg/s, 
bright AGN: Lagn > 10^^ erg/s 


green valley, in the distribution of AGN. A detailed compar¬ 
ison between active and inactive galaxies could potentially 
reveal differences in the two galaxy populations, but this 
exercise is beyond the scope of this study. We finally note 
that the Matar — SFR plane at higher redshifts, {z ~ 2) is 
characterised by the absence of the passive regime (for both 
active and inactive galaxies). The bimodal distribution of 
AGN at 2 = 0 is linked to the quenching of galaxies in the 
low-redshift universe and relates strongly to the colour bi¬ 
modality of AGN in GALFORM, as shown in|Georgakakis et al.| 
|M4|. 


To gain more insight into the SF properties of AGN 
hosts, we show in the top panels of Fig. the molecular 
hydrogen gas content as a function of AGN luminosity for 
a small subset of galaxies at 2 = 0 and 2 = 2.1. We remind 
the reader that the molecular hydrogen content is calculated 


following Lagos et al. (20111, where the molecular-to-total 


gas ratio increases with the mid-plane hydrostatic pressure 
of the galactic disk. What is immediately evident in this 
plot is that both at high and low redshifts the hosts of faint 
and bright AGN reach similarly high molecular gas contents, 
but with a larger scatter towards lower values in faint AGN 
hosts at 2 = 0. At a first look, this is surprising given that the 
hosts of faint AGN tend to be more passive systems. Indeed, 
the total cold gas reservoirs of early-type galaxies are lower 
than those of starbursts in our model. Yet, the model pre¬ 
dicts a strong correlation between the molecular-to-atomic 
hydrogen ratio with increasing bulge-to-total ratio, meaning 
that early type galaxies are relatively richer in molecular 
hydrogen compared to atomic hydrogen, as implied also by 


past and recent observations (jYonng & Knezek 

1989 

Bettoni 

et al. |2003| |Leroy et al.||20081 |Lisenfeld et al. 

[ 2 M 1 

IBoselli 


et al. 20141. This is due to the higher compactness of early- 


type galaxies compared to starburst galaxies, which results 
in higher gas pressure and thus, higher molecular-to-atomic 
hydrogen ratios ( Lagos et al.|2014b (. 

Despite the high H 2 mass in early-type galaxies, the 
SFR in these systems remains low compared to starburst 


galaxies, as has already been shown in the lower panel of 
Fig-i At high redshift this is due to the different timescales 
for star formation between spheroidals and starburst sys¬ 
tems, with the latter ones having shorter timescales (see also 
section 2.31. At low redshifts, early-type galaxies, i.e. bulge- 


dominated systems, tend to have relatively low mean SFRs 
due to their lower mean molecular gas contents. This is also 
shown in the lower panel of Fig. where we now plot the 
specific SFR (sSFR). 

Similarly to the SFR, the AGN luminosity is also low 
in early-type galaxies. This is due to the fact that these sys¬ 
tems are dynamically stable, having very subdominant disks. 
They rarely undergo disk instabilities and thus, quasar-like 
AGN activity is never achieved (though major mergers could 
occasionally produce significant accretion). BHs typically 
grow via the hot-halo mode, which in the redshift regime 
of interest always produces low density accretion flows, irre¬ 
spective of the cold gas reservoir of the galaxy. In contrast, in 
the high-luminosity regime we find gas-rich starburst galax¬ 
ies that often experience disk instabilities which then trigger 
efficient and quasar-like accretion onto their central BHs. 

Our results suggest that galaxies in the low luminosity 
region of the Lagn — fLeonm plane are less efficient both at 
building stars and growing BHs than their high luminosity 
counterparts. 


4 DISCUSSION 


Overall our model produces a SFR-LAGNrelation that is in 
reasonable agreement with observations, especially taking 
into account that no parameters of the model were tuned 
to reproduce this specific observational relation. As it can 
been seen in Fig. the model predicts FIR luminosities for 
the faint and moderate luminosity AGN that match those 
observed in the Herschel surveys. In addition, the model pro¬ 
vides a reasonably good fit to the FIR luminosity function of 
the total galaxy (Lacey et al. in prep) and AGN population 
( [Fanidakis et al.]|2012[ ). Note that [Rosario et ^ use stack¬ 
ing to estimate the FIR flux for sources whose flux is below 
the 3 (t detection limit of the RAGS instrument. Therefore, 
there is in principle no reason to create a mock catalogue for 
mimicking flux limits and biases related to the observations. 
Thus, assuming that their data sample is complete, possible 
reasons for the disagreement could be the following. First 
of all, there could be a strong contribution in the model of 
very faint FIR sources to the average value. 

As we have already mentioned, the objects populating 
this region of the plot are identified with early type galax¬ 
ies. These systems are passive in terms of their efficiency in 
forming of stars (although occasionally a starburst could oc¬ 
cur as a consequence of a major merger), and therefore their 
SFR and, as a consequence, their FIR emission depend on 
the H 2 gas reservoir ( Lagos et al.|20lT |. However, the model 
seems to produce a few too many low mass systems that are 
also low in cold gas (atomic and molecular) and thus have 


low SFR. But based on previous work with this model (La- 
gos et al.|2012 2014a|b[ also Lacey et al. in prep.), the FIR 
luminosities and SFRs in the vertical axes of Fig. [^should 
be in principle consistent with the observations. 

A more plausible reason for the strong disagreement 
could be the simplicity of the assumptions for the triggering 
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Figure 6. Top: Predicted molecular hydrogen mass, Mjj^, as a function of bolometric AGN luminosity for 4,000 randomly selected 
galaxies at 2: = 0 (left) and 2 = 2.1 (right). Galaxies in which the central BH experiences starburst-mode activity are shown in blue, while 
those with BHs growing during the hot-halo mode in red. Bottom: Predicted specific SFR (SFR/M*) as a function of AGN luminosity. 
Golors and redshifts as at the top panel. We only show galaxies with stellar masses above lO^'^ M©. 


and variability of AGN activity. In GALFORM, AGN activ¬ 
ity begins with the onset of star formation in a starburst 
without a time delay. A possible delay between the trigger¬ 
ing of AGN activity and the formation of stars in a burst 
could result into a decrease towards lower AGN luminosities 
for the objects populating the bright AGN regime of the 
diagram. The incorporation of an AGN variability model 
could have a similar effect. AGN luminosities are calculated 
assuming a constant accretion rate. Thus, the luminosity 
is constant during the entire course of accretion. However, 
AGN are known to exhibit strong variability on a wide range 
of timescales ( [Novak et al.|20fT| . 

Recently, Hickox et al. presented a model in which SF 
and BH accretion are closely connected over long timescales, 
but that short-term AGN variability can wash out this cor¬ 
relation for low to moderate Tagn- Here, despite the fact 
that the GALFORM intrinsic FIR-AGN relation is driven by 
completely different phenomena, the inclusion of AGN vari¬ 
ability could improve the agreement with the observations 
simply by shifting some bright quasars to lower AGN lumi¬ 
nosities. 

A final explanation for the disagreement with the ob¬ 
servations at lower AGN luminosity could be that the AGN 
luminosities we calculate for the AGN in the hot-halo mode 
are underpredicted (see also discussion about the hot-halo 


AGN in Krumpe et al. in prep). Indeed, the sharp transi¬ 
tion from the hot-halo mode to the starburst mode, which 
appears as a strong break at Lagn ~ 10^^ erg s“^ in the 
Tagn — correlation, could possibly be smoothed out 

if the hot-halo AGN luminosities were systematically higher. 
In the current version of the model, the hot-halo luminos¬ 
ity is calculated directly from the cooling properties of the 
host DM halo, via the expression M = I/cooi/(ekinC^). While 
Lcooi is well defined in the model, the value of the efficiency 
parameter Ckin is loosely constrained, mainly by requiring 
the model to reproduce the BH mass function in the local 
Universe. By boosting the luminosity of hot-halo AGN we 
could in principle make the transition from one regime to the 
other smoother and thus obtain a better agreement with the 
data. We note however that our aim in this paper is merely 
to report what observable correlations are predicted by a 
model in which BH and SFR are strongly coupled. We made 
no adjustments to the model in preparation for doing this 
analysis. Nevertheless, it is interesting that a comparison 
with FIR observations of X-ray selected AGN could provide 
possible constraints in the modelling of AGN feedback in 
galaxy formation models. 
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5 CONCLUSIONS 


In the current paradigm of galaxy formation the triggering 
mechanisms of AGN activity are often responsible for ini¬ 
tiating intense SF in the host galaxy. At the same time, 
it is widely accepted, yet unproven, that AGN activity is 
the main driver of the quenching of SF in massive galax¬ 
ies (Di Matteo et al. 2005| Monaco & Fontanot 20051. In 
both cases, a correlation between observable proxies for SF 
and AGN activity is somehow expected. However, several 
recent observational studies have shown that there is little 
I Rosario et al.|2012 1 or practically no correlation (Mullaney 


et al. 2012b Stanley et al. 20151 between FIR luminosity 


and the AGN luminosity. These results have been seen as a 
challenge for current galaxy formation models. 

In this study, we have reported the predictions of the 
semi-analytic model GALFORM for the connection between 
SFR and AGN activity. GALFORM calculates AGN proper¬ 
ties by tracking the BH accretion rate and the evolution of 
BH mass and spin. BHs grow via cold gas accretion (star- 
burst mode), usually triggered by disk instabilities or galaxy 
mergers, and by hot-gas accretion (hot-halo mode), typically 
in massive quasi-hydrostatic haloes that are subject to AGN 
feedback. During the course of accretion the GALFORM code 
calculates various AGN properties, as for example the disk 
emission in different bands, by considering the dependence 
of the bolometric disk emission on the structure of the ac¬ 
cretion disk. 


Galaxies in GALFORM build stars in disks from molecular 
hydrogen. When a galaxy experiences a disk instability or 
galaxy merger all the available cold gas (atomic and molec¬ 
ular) turns into stars via a burst. The growth of BHs is 
strongly linked to the buildup of stars. During a burst of SF 
a fraction of the gas that turns into stars is accreted onto the 
BH (starburst mode). This creates a strong correlation be¬ 
tween SFR and AGN luminosity. In galaxies found in more 
massive, quasi-hydrostatic haloes SF is less efficient. The ac¬ 
cretion power couples with the cooling properties of the gas 
in the halo resulting in a suppression of gas cooling and a 
regulation of the SF. This gives raise to an anti-correlation 
between SFR and AGN activity. 

To compare our predictions to recent Herschel observa¬ 
tions for the mean SFR of X-ray selected AGN we compute 
the properties of galaxies in the FIR. The FIR emission in 
our model is due to the reprocessing by dust of the inci¬ 
dent stellar continuum radiation. The resulting emission at 
60 pm, uLsoiim, scales nearly linearly with SFR (Fig.[^, and 
thus represents a good tracer for the SFR. For the brightest 
AGN, GALFORM predicts a strong correlation between Lboi 
and vLeoum, which arises from the strong coupling of the 
BH accretion to starbursts. For faint and intermediate lu¬ 
minosity AGN, the model predicts a mildly negative correla¬ 
tion, which is a consequence of the negative feedback effect 
these AGN have on their host galaxies. When compared to 
the Herschel PACS data for X-ray selected AGN in the COS¬ 
MOS and GOODS-S/N fields we find a very good agreement 
in the bright AGN luminosity regime. In the low-luminosity 
regime we find that the model systematically underpredicts 
the average u Leo lira- 

Finally, we showed that the objects populating the 
bright and faint regimes of the Lboi — vLeoiim plane represent 
different classes of galaxies. In the high luminosity regime 


we find gas-rich disk galaxies, that recently underwent a disk 
instability, and now exhibit strong starbursts and prodigious 
BH growth and AGN activity. In the low luminosity regime, 
we find massive early-type galaxies that experience quiescent 
SF, BH growth and AGN activity. Interestingly, despite the 
fact that these galaxies are subject to AGN feedback, many 
reach relatively large masses of molecular hydrogen. In fact, 
we find that their H 2 reservoirs at high redshift are similar 
to the gas-rich disk galaxies in the high luminosity regime 
of the diagram. These galaxies, however, are inefficient in 
building stars as can be seen in their sSFRs. We find that 
their SFRs and FIR luminosities tend to be on average ap¬ 
proximately one order of magnitude lower than the galaxies 
in the bright luminosity regime. Nevertheless, a large frac¬ 
tion reach high molecular hydrogen masses, so a large popu¬ 
lation of low AGN luminosity galaxies should have as bright 
CO emission as their high-AGN luminosity counterparts. 

Future observations of the molecular gas/CO emission 
in X-ray selected AGN will help to further disentangle the 
SFR-AGN correlation. Until then, we hope that this work 
will provide adequate motivation that AGN activity remains 
a viable solution to explain the origin of red-and-dead galax¬ 
ies. 
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